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FOREWORD 


This  report  describes  a  portion  of  the  results  obtained  on 
NASA  Grant  NSG  3044.  This  work  was  done  under  subcontract  to  the 
University  of  Illinois,  Urbana,  with  Prof.  S.S.  Wang  as  the  Prin¬ 
cipal  Investigator.  The  prime  grantee  was  the  Massachusetts  Insti¬ 
tute  of  Technology,  with  Prof.  F.J.  McGarry  as  the  Principal  Inves¬ 
tigator  and  Dr.  J.F.  Mandell  as  a  major  participant.  The  NASA  - 
LeRC  Project  Manager  was  Dr.  C.C.  Chamis. 

Efforts  in  this  project  are  primarily  directed  towards  the  de¬ 
velopment  of  finite  element  analyses  for  the  study  of  flaw  growth  and 
fracture  of  fiber  composites.  The  analysis  of  such  problems  using 
three-dimensional  analyses  may  be  limited  by  the  presence  of  very 
localized,  high  stress  gradients  as  occur  at  free  edges.  The  work 
described  in  this  report  is  a  theoretical  investigation  of  such  ef¬ 
fects  resulting  from  thermal  or  hygroscopic  loading.  The  results 
given  here  are  for  thermal  effects,  but  the  same  treatment  also  ap¬ 
plies  to  hygroscopic  effects,  differing  only  by  the  expansion  co¬ 
efficient  used.  Analogous  results  to  those  given  here,  but  for  hy¬ 
groscopic  loading,  may  be  found  in  the  following  papers: 

(1)  S.S.  Wang  and  I.  Choi,  AIAA  paper  80-0713-CP,  21st  Proceedings 
of  AIAA/ASME/ASCE/AHS  Structures,  Structural  Dynamics  and  Mate¬ 
rials  Conf.,  Seattle,  WA,  1980;  also  to  appear  in  AIAA  Journal, 
November,  1982. 

(2)  S.S.  Wang  and  I.  Choi,  in  Developments  in  Composite  Materials 


and  Structures,  J.R.  Vinson,  ed.,  ASME ,  p.  315  (1980). 


ABSTRACT 


Thermal  stresses  near  geometric  boundaries  of  fiber-reinforced  composite 
laminates  play  an  Important  role  in  controlling  complex  failure  modes  and 
ultimate  performance  of  composite  materials  in  severe  thermal  environment. 

This  paper  presents  an  investigation  of  boundary-layer  thermal  stress  singularity 
and  distributions  in  composite  laminates.  Based  on  the  theory  of  anisotropic 
thermoelasticity,  a  system  of  coupled  governing  partial  differential  equations 
is  obtained.  Edge  boundary  conditions  and  interface  continuity  conditions  lead 
to  a  transcendental  characteristic  equation  for  determining  thermal  stress 
singularity  at  the  laminate  boundary.  Complete  thermal  boundary-layer  stress 
and  displacement  solutions  are  obtained  by  an  eigenfunction  expansion  method 
in  conjunction  with  a  boundary  collocation  procedure.  The  thermoelasticity 
solution  in  the  region  away  from  the  singular  domain  is  found  in  excellent 
agreement  with  existing  approximate  numerical  results.  As  the  edge  is  approached, 
the  singular  terms  control  the  near-field  behavior  of  thermal  boundary-layer 
stresses.  Results  are  presented  for  the  cases  of  various  angle-ply  graphite-epoxy 
laminates  with  [9/-9/-0/0]  fiber  orientations.  Thermal  boundary-layer  thickness 
is  defined  by  considering  strain  energy  density  distribution  along  ply  interface. 
The  thermal  boundary-layer  thickness  is  shown  to  depend  on  the  anisotropy  of 
individual  lamina,  ply  thermomechanical  properties,  and  relative  thickness  of 


adjacent  layers. 
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1.  INTRODUCTION 

The  response  of  a  composite  laminate  near  its  geometric  boundaries  sub¬ 
jected  to  severe  thermal  and  other  environmental  loading  has  attracted  much 
attention  recently,  since  failure  of  composite  materials  is  frequently  initiated 
at  the  boundaries.  The  thermal  stress  field  in  the  vicinity  of  the  laminate 
boundaries,  which  may  be  primarily  responsible  for  strength  degradation  and 
failure  of  composites,  has  been  Investigated  by  several  researchers  using 
different  approximate  methods  [1-4].  Approximate  solutions  for  the  thermal 
boundary-layer  problem  have  revealed  several  unusual  features.  Interlaminar 
stresses  near  a  traction-free  boundary  of  a  composite  laminate  subjected  to  a 
uniform  temperature  change  have  been  found  to  be  very  high  and  inherently 
three-dimensional.  It  has  also  been  reported  that  the  high  thermal  stress  is 
confined  within  a  localized  region  of  several  laminar  thicknesses  from  the 
edge  and  that  in  the  boundary-layer  region  they  cannot  be  assessed  accurately 
by  the  classical  lamination  theory  [5,6].  The  behavior  of  this  highly  stressed 
boundary-layer  region  is  of  great  Importance  in  controlling  complex  failure 
modes  and  ultimate  performance  of  the  composites.  Understanding  the  fundamental 
nature  of  boundary-layer  thermal  stresses  is  essential  to  the  failure  analysis, 
design,  and  processing  optimization  of  composite  materials. 

While  all  previous  numerical  approximate  solutions  indicated  an  unbounded 
trend  of  boundary-layer  thermal  stresses  and  postulated  possible  existence  of 
a  stress  singularity  at  the  edge,  the  search  for  the  order  or  strength  of  the 
stress  singularity  has  been  unsuccessful.  No  information  concerning  the  exact 
order  of  the  boundary-layer  stress  singularity  has  been  reported  yet,  to  the 
authors’  knowledge.  The  apparent  difficulty  may  result  from  the  complicated 
nature  of  the  problem  such  as  the  therraomechanical  anisotropy  of  each  individual 
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fiber-reinforced  lamina,  the  geometric  discontinuity,  and  the  abrupt  change  of 
material  properties  through  the  laminate  thickness  direction.  Since  the 
thermal  boundary-layer  effect  is  localized  in  nature,  it  is  necessary  to 
determine  the  exact  order  of  the  laminate  edge  stress  singularity  so  that 
complex  thermal  response  in  the  vicinity  of  laminate  boundaries  can  be  studied 
more  accurately.  This  paper  presents  a  rigorous  theoretical  study  of  thermal 
boundary-layer  stress  singularity  and  distributions  in  composite  laminates 
subjected  to  uniform  thermal  loading. 
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2 .  FOEMULATION 

2 . 1  Basic  Equations 

Consider  a  general  composite  laminate  composed  of  fiber-reinforced 
plies  subjected  to  mechanical  and  thermal  loading,  as  shown  in  Fig.  1.  Denote 
the  constitutive  equation  of  each  individual  ply  by  the  Duhamel-Neumann  form 
of  generalized  Hooke's  law  in  contracted  notation  as 


e.  =  S..  a.  +a.  AT, 
1  13  3  1 


(i,j  =  1,2,3, ... .6) , 


(1) 


where  the  repeated  subscript  indicates  summation;  is  the  compliance  tensor; 
a^,  the  thermal  expansion  coefficient,  and  AT,  the  change  of  temperature.  The 
engineering  strains,  e,,  in  Eq.  1  are  defined  in  a  Cartesian  coordinate  system 
by 
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where  u,  v  and  w  are  displacement  components.  The  stresses,  are 
defined  in  an  analogous  manner. 

The  composite  laminate  considered  here  has  a  finite  width  and  is  sub¬ 
jected  to  a  uniform  axial  extension,  e,  along  the  z-axis  and  a  uniform  temper¬ 
ature  change,  AT.  The  composite  is  assumed  to  be  sufficiently  long  that,  in 
the  region  far  from  the  end ,  the  end  effect  is  neglected  by  virtue  of  Saint 
Venant's  principle.  Consequently,  stresses  in  the  laminate  are  independent 
of  the  z  coordinate.  The  special  case  where  e  vanishes  identically  corresponds 
to  the  well  known  generalized  plane  deformation  [7].  Under  these  assumptions, 
equations  of  equilibrium  without  body  force  read 
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(3c) 


Integrating  Eq.  1  with  the  aid  of  Eq.  2,  one  can  obtain  u,  v  and  w  as 


2  5  9W 

u  .  -  (Sj  “5  “  -  -ST*  VjK.y), 


(4a) 


~2  ;\T) 

V  =  -  ^  |-+  (S.  .  a.  +  a,  AT 
2  9y  4j  j  4 
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17)  ^ 


(4b) 


w  =  D  z  +  W^(x,y)  , 


(4c) 


where 


D  =  S„.  a.  +  a-  AT, 
3j  J  3  ’ 


(5) 


and  U  ,  V  and  W  are  arbitrary  functions  of  x  and  y  only. 

000 

Following  the  procedure  in  [7],  it  can  be  shown  after  some  mathematical 
manipulation  that  general  expressions  for  the  displacements  and  the  stress 
component  have  the  following  forms : 


u  =  -  J  yz  +  U(x,y)  +10^  z  -  y  +  u^,  (6a) 

V  =  -  Y  ^2  ^33  xz  +  V(x,y)  +  x  -  z  +  v^,  (6b) 

w  =  (A^  X  +  y  +  A^)S22  z  +  W(x,y)  +  co^  y  -  0)2  x  +  w^,  (6c) 

X  +  A2  y  +  A^  ~^^3j  ^  “s  ^3  =  1,2, 4, 5, 6).  (6d) 

The  unknown  functions,  U,  V  and  W,  depend  on  x  and  y  only,  and  can  be  shoxm 
easily  to  obey  the  following  relationships: 


9x 


Oj  +  S^3(A^  X  +  A2  y  +  A^)  +  AT, 


(7a) 


5 


=  S2J  +  S23(Aj^  X  +  A2  y  +  A^)  +  AT,  (7b) 

3W  ~ 

^  =  Sj  ‘"j  h3^h  X  +  A2  y  +  A3)  +  y  +  AT.  (7c) 

3W  ~ 

"  ^4j  ""j  ^  +  ^2  y  +  A3)  -  X  +  AT,  (7d) 

9^  +  ^  =  ^6j  X  +  A2  y  +  A3)  +  AT,  (7e) 


where 


(j  =  1,2, 4, 5, 6). 


S.  . 
ij 


S.  . 

13 


~  ^13  ^j3'^^33’ 


(8a) 


“i  "  “i  “  “3  ^13^^33  1,2, 4, 5, 6).  (8b) 

It  is  obvious  that  the  constants,  u^,  v^,  and  w^(i  =  1,2,3),  in  Eqs.  6a-6d 
characterize  the  rigid  body  translation  and  rotation  of  the  solid.  A^^  and 
represent  the  bending  of  the  laminate  in  the  x-z  and  y-z  planes,  respectively. 
A3  characterizes  the  uniform  axial  extension  of  the  composite  laminate,  and  A^, 
the  relative  angle  of  rotation  about  the  z-axis . 


2.2  Govemi-nff  Partial  DifferentiaZ  Equations ^  and  Boundary  and  End  Conditions 
Introducing  Lekhnitskii’s  stress  functions,  F  and  'P,  such  that 


a  = 


dJ_F 
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32  F 
3x2  > 


xy 


3^  F 
3x  3y’ 


xz 
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T  =  - 

yz 


3y 

3x’ 


(9a-e) 


one  can  show  that  the  equations  of  equilibrium  are  satisfied  identically. 
Eliminating  U  and  V  from  Eqs.  7a,  7b  and  7e,  and  W  from  Eqs.  7c  and  7d,  one 
obtains  the  following  system  of  governing  partial  differential  equations: 
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F  +  L2  W  2A^  +  A^  AT, 


F  +  -  H2  AT, 


(10a) 

(10b) 


where  L^,  L.,  L, ,  H, ,  and  H„  are  linear  differential  operators  defined  as 
2  3  4  1  2 
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Now  consider  boundary  conditions  on  the  lateral  surfaces  and  at  the  ends 
of  the  laminate  strip.  Assuming  that  the  lateral  surfaces  are  free  from  sur¬ 
face  tractions,  one  may  have  the  following  conditions: 


and 


a  n  +  T  n  =  0, 
XX  xy  y 


T  n  +  a  n  =0, 
xy  X  y  y 


T  n  +  T  n  =  0, 
xz  X  yz  y 


(12a) 

(12b) 

(12c) 


where  n  and  n  denote  the  directional  cosines  of  the  outward  unit  normal  on 
X  y 

3B  as  shown  in  Fig.  2.  The  conditions  at  the  ends  of  the  composite  laminate 
be  obtained  from  the  statically  equivalent  loading  as 


may 
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T  dx  dy  =  0 
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JJ 
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T  y)  dx  dy  =  M^, 
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(13a) 


(13b) 


(13c) 


(13d) 


(13e) 


(13f) 


at  z  =  ±L,  where  P^,  M^,  are  the  applied  axial  force,  bending  moments, 

and  twisting  moment,  respectively. 


2.3  Interface  Continuity  Conditions 

Consider  a  portion  of  the  laminate  composed  of  two  different  fiber- 
reinforced  laminae,  as  shown  in  Fig.  2.  Stress  functions  of  Eqs.  9a  and  9b 
must  be  defined  for  each  ply.  Assuming  that  the  plies  are  perfectly  bonded 
along  the  interface  BB^.,  one  can  immediately  establish  the  continuity  con¬ 
ditions  of  the  stress  and  displacement  along  the  interface  between  the  kth  and 
(k+l)th  plies  as  the  following: 


,(k) 

+ 

P 

II 

_^(k+l) 

M 

1 

„(k+l) 

(14a) 

X 

X 

xy 

y 

X 

X 

xy 

y 

(k) 

(k) 

n 

+ 

_  ^(k+1) 

„(k+l) 

- 

„(k+l) 

(14b) 

xy 

X 

y 

y 

xy 

X 

y 

y 

(k) 

(k) 

n 

+ 

x»> 

n(^)  = 

_ 

_^(k+l) 

_  x^k+i) 

(k+l) 
n  , 

(14c) 

xz 

X 

yz 

y 

xz 

X 

yz 

y 

8 


^(k)  ^  ^(k+1)^ 
^(k)  ^  ^(k+1)^ 


(15a) 

(15b) 

(15c) 


By  using  divergence  theorem,  it  can  be  readily  seen  that  Eqs.  13a  and  13b  are 
satisfied  identically  by  virtue  of  3c,  12c  and  14c.  Substituting  the  displace¬ 
ments,  Eqs.  6a-c,  into  the  above  continuity  conditions,  Eqs.  15a-c,  one  finds 
that  the  unknown  constants  for  two  adjacent  layers  in  Eq.  6  are  related  by 


o  o 
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(16a-c) 
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and 
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3.  METHOD  OF  SOLUTION 

The  complex  governing  differential  equations  and  boundary  conditions 
formulated  in  the  previous  section  may  be  simplified  considerably  as 
appropriate  loading  conditions  and  geometric  symmetry  are  taken  into  con-  . 
sideration.  For  simplicity  and  without  loss  of  generality,  the  following 
geometric  and  environmental  conditions  are  introduced  for  the  present  thermal 
boundary-layer  stress  problem: 

(a)  the  temperature  change  AT  is  constant  and  uniformly  distributed 
throughout  the  laminate,  and  there  is  no  external  mechanical  loading  applied; 

(b)  the  interface  is  a  straight  line  and  meets  the  traction-free  edge  of 
the  composite  laminate  by  a  right  angle. 

Under  these  assumptions  the  governing  differential  Eqs.lOa  and  10b  may 
be  simplified  as 

F  +  L2  '1'  =  -  2A^  +  -  A2  (17a) 

< 

F  +  f  =0,  (17b) 

with  boundary  and  end  conditions, 

a=-r=T=0,  onx  =  0,(18a-c) 
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on  B, 


(19a) 


(x  X  -  T  y)dx  dy  =  0,  on  B.  (19b) 

^  yz  xz-^ 

B 

The  Interface  continuity  conditions  of  tractions  and  displacements  may  be  ex¬ 
pressed  as  the  following: 
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„(k)  __  p(k+l) 


.  -  V 


(k+l)^(k+l) 


„(k)  .  „(k+l)  (20b) 


along  y  =  0,  as  the  original  of  the  coordinate  system  is  moved  to  the  free  edge. 


3.1  '  Eomogeneous  Solution 

The  simplified  governing  differential  equations,  Eqs.  17a  and  17b,  are 
coupled,  linear  partial  differential  equations  with  constant  coefficients  re¬ 
lated  to  the  anisotropic  thermomechanical  elastic  constants  of  each  composite 
lamina.  The  general  solution  consists  of  a  homogeneous  solution  and  a  par¬ 
ticular  solution  depending  on  geometric,  loading,  and  boundary  conditions  of 
the  problem  under  consideration.  Lehknitskii  has  shown  [7]  that  the  homo¬ 
geneous  solution  of  the  above  mentioned  system  of  governing  partial  differential 
equations  has  the  general  form 


6 

F(x,y)  =  y  F,  (x  +  p  y),  (21a) 

k=l 


6 

’i'(x,y)  =  I  \  \ 

k=l 


(21b) 


where  the  prime  (')  in  Eq.  21b  denotes  differentiation  of  the  function 
Fj^(x  +  y)  with  respect  to  its  argument,  and  the  coefficients  Pj^  are  the 
roots  of  the  following  algebraic  characteristic  equation: 


and 


with 


^^(y)  ^2^^^  “  ^3(y)  =  0, 


^2^1^)  S55  p^  ^^45  ^  ^44’ 


(22a) 

(22b) 

(22c) 
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Jl3(y)  =  S33  p3  _(s^^  +  S3g)p2  +(833  +  p  -  83^,  (22d) 

^^(P)  =  S33  p**  -  28^g  p3  +(28^2  +  8gg)  p2  -  282^  P  +  822-  (22e) 


It  has  been  shown  [7]  that  Eq.  22a  cannot  have  a  real  root.  Thus  the  roots  p^^ 
appear  as  complex  conjugates,  and  Fj^  are  analytic  functions  of  the  complex 
variables  Zj^  =  x  +  Pj^  y.  8ubstituting  the  expressions  of  F(x,y)  and  ’i'(x,y) 
in  Eqs.  21a  and  21b  into  Eqs.  9a-e,  the  homogeneous  components  of  the  stress 
a.  may  be  expressed  in  terms  of  F,  (Z,  )  as 

1  iC  1C 


“f’  ■  J,  “k 

k=l 


y  k=i 


(h)  =  _  6 


(h)  _ 


•‘k  ’’k  ''k<V> 

k=l 


,(h) 

xy 


k=l 


(23a) 

(23b) 

(23c) 

(23d) 

(23e) 


The  expressions  for  displacement  components  may  be  obtained  directly  from 
Eqs.  7,  10  and  23  with  omission  of  terms  which  are  to  be  included  in 
the  particular  solution.  Hence,  one  has 


i  ■’k  W- 


k=l 


■  I  %  W- 


k=l 


(24a) 


(24b) 
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and 

=  J,  W  ■ 

k=l 

where 

Pk  =  ^11  Pk  +  ^12  -  ^14  Pk  +  ^15  Pk  Pk  -  ^16  Pk’ 


\  =  ^12  Pk  +  ^22/Pk  ■  ^24  Pk^Pk  +  ^25  Pk  -  ^26’ 


\  ^14  Pk  ^24  ^Pk  ^44  Pk'^Pk  ■*■  ^45  ^k  "  ^46  ‘ 


We  now  choose  the  form  of  F,  (Z,  )  as 

k  k 


\(\)  =  %  Z^'^^/[(6+l)(6+2)], 


where  and  6  are  arbitrary  complex  constants  to  be  determined  later, 
stituting  Eq.  25  into  Eqs.  23  and  24  gives 


,,(h)  _  Y  ro  ,,2  7'^  4.  r  TT^  7^1 

^  Pk  \  S+3  Pk 

k=l 


oi*'’ -  + 


k=l 


,(h) 

yz 


-  "k<+V3 

k=l 


and 


=  y  [C,  n,  y,  zf  +  c,  n,  y, 
xz  ,  ^  T  k  k  k  k  k+3  k  k  k 
k=l 


^(h) 

xy 


^  Pk  \  \+3  Pk 


k=l 


t"k  Pk  +  "k4-3  Pk 

k=l 


(24c) 

(24d) 

(24e) 

(24f) 

(25) 

Sub- 

(26a) 

(26b) 

(26c) 

(26d) 

(26e) 


(27a) 
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-  J,  '‘^k  %  4'"'  +  V3  \  C' 

k=l 

*  ,  i  '‘^k  '=k  C  ^  ^3  \  .  (2^<=) 


where  the  overbar  denotes  the  complex  conjugate  of  the  associated  variable. 

The  homogeneous  solutions  for  the  stress  and  displacement  shown  in 
Eqs.  26  and  27  are  required  to  satisfy  homogeneous  boundary  conditions 
and  Interface  continuity  conditions,  i.e.,  Eqs.  18  and  20.  This  leads  to  a 

standard  eigenvalue  problem,  which  results  in  a  complicated  transcendental 

equation  for  determining  the  value  of  6.  By  taking  the  complex  conjugate 

of  Eqs.  26  and  27,  it  is  readily  seen  that  always  appear  as  pairs  of 

complex  conjugates.  Thus  the  stresses  and  displacements  can  be  made  real  by 

superposing  the  conjugate  solutions.  It  is  noted  that  a  similar  method  of 

extracting  the  real  function  for  an  isotropic  material  case  was  discussed  by 

Theocaris  [8]  who  assumed  the  stress  function  including  both  6  and  6  in  the 

formulation.  It  is  also  noted  that,  due  to  positive  definiteness  of 

strain  energy  of  an  elastic  body,  the  value  of  6  bounded  by  -1  <  Re [6  ]  <0 

m  m 

characterizes  the  order  of  singularity  of  the  boundary-layer  stress  field  in 
the  composite  laminate.  It  has  been  found  in  a  related  paper  [9]  that  there 
is  only  one  real  6^  which  meets  this  condition  in  the  singular  boundary-layer 
stress  problem.  Detailed  information  on  the  edge-stress  singularity  has  been 
given  in  Reference  9. 
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(j  =  1.2,3), 


(28b) 


where  the  unknowns  d  are  real,  and  f.  and  g.  denote  known  eigenfunctions 

m  im  jm 

corresponding  to  the  m-th  eigenvalue  6^. 

3.2  Farticutar  Solution 

The  particular  solution  for  the  problem  may  be  sought  in  the  form  of 
polynomials  as 

(x,y)  =  a^  +  a2  y  +  xy^  +  a^  y^  +  a^  x^  +  a^  xy  +  a^  y^ ,  (29a) 

(x,y)  =  ag  x2  +ag  xy  +  a^^  y2  +  a^^^^  x  +  a^2  y>  (29b) 

where  a^  are  arbitrary  constants  to  be  determined.  It  is  seen  that  Eq.  17b 
is  satisfied  identically  and  that  Eq.  17a  provides  the  following  relation  for 
each  ply: 


^^24  ^1  ^^^25  ^46^^2  “  ^^^14 

+  63^3  a^  +  23^^  ag  -  23^^  a^  +  23^^  a^^ 

=  -2A^  +  3g^  -  A2  3^3. 

(30) 

3ubstituting  Eqs.  29a  and  29b  into  Eqs.  7  and  9  gives 

=  2a„  X  +  6a,  y  +  2a,,, 

X  3  4  7 

(31a) 

=  6a,  X  +  2a„  y  +  23^, 
y  1  2  5 

(31b) 

(p) 

-  -  2ag  X  -  a,  y  -  a^^. 

(31c) 

"xa’  -  ^  +  2a^0  ='  +  ‘l2- 

(31d) 

-  232  X  -  2a3  y  -  =6’ 

(31e) 
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3U 


(P) 


3x  .  ^11  +  ^13’ 


(32a) 


3V 


(p) 


3y  '  ^21  ^  ^22  ^  '*'  ^23’ 


(32b) 


3W 


(P) 


3W 


3y 

(P) 


=  (E^^  -  A^)  X  +  E^2  y  +  ^43’ 


3x  =  ^51  ^  +  ^^2  +  \>  y  +  ^53’ 


3U^P^  ^  3V^P^ 


3y 


3x  ■  =  ^1  ^  +  ^62  y  +  ^63’ 


where 

E  _  =  2S  _ 

+ 

6S.„ 

a. 

-  2S,,  a„  +  S 

jl 

jl 

3 

j2 

1 

j4  8 

E  „ 

=  6S  , 

a,  + 

2S  „ 

a^ 

-  S.,  a^  +  2S 

j2 

jl 

4 

j2 

2 

j4  9 

E  ^ 

=  2S., 

a_  + 

2S.^ 

—  S.,  a, -  +S 

j3 

jl 

7 

j2 

5 

j'^i  11 

(32c) 

(32d) 

(32e) 

(33a) 

(33b) 


Integrating  Eqs.  32a-e  with  the  aid  of  Eq.  30,  one  can  obtain 


U 


2  ®11  +^12^'*'  ^13  ^  2^^62  ~  ^21^  y^  2  ^63 


(34a) 


V  ^  -  ^21  2  ^22  y^  ^  y  ^  ^  "^9 


23 


W 


2'  61  12^ 
1 


2  63 


j  x2  +  (E^2  +  A^)  xy  +  E^2  X  +  I  E^^  +  ^43  y- 


(34b) 

(34c) 


Thus  the  particular  solution  for  the  displacement  can  be  written  as 


u 


— A  S  z^-A 
2  ^1  ^33  ^  4 


.(P)  -  1 


2  2  33 


^^\x,y)  +  0)^  z  -  u)^  y  +  u^. 

(35  a) 

*‘^^(x,y)  +  X  -  z  +  v^. 

(35b) 

w^P^=  (A^  X  +  A^  y  +  A^)  z  +  W^^^x.y)  +  y  -  x  +  w^ . 


(35c) 
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Equations  31  and  35  are  required  to  satisfy  the  boundary  conditions,  Eqs .  18a-c, 
and  the  interface  continuity  conditions,  Eqs.  20a-b  of  the  current  problem. 

This  leads  to  the  establishment  of  the  following  relations 


(m) 

=  a}“^  = 

=  a^“^  = 

aj“^  = 

=  0 

3 

4 

6 

7 

10 

12 

(k) 

-  Jk+1) 

„(k+l) 

-  Jk+1) 

-  , 

^2 

-  ^2 

’  ^5 

^5 

.(k) 

^8 

-  ^8 

,  ^11 

-  hi 

(k) 

E(k) 

-  E^’^)  = 

11 

11 

13 

13 

’  61 

12 

(36b-f) 


(k+1)  (k+1),  . 

^61  "  ^12 


1  (k)  (k)  _  1  p(k+l)  (k+1)  (k)  _  (k+1)  (k)  _  (k+1) 

2  ^63  +  “3  “2  ^63  ^  ^3  ’  hi  ~  ^51  ’  ^53  “  ^3  *  (36j-l) 


By  examining  Eqs.  16,  30  and  36,  it  is  observed  that  there  are  44  unknowns  (in¬ 
cluding  the  unknown  rigid-body  translations  and  rotations)  related  by  34  linear 

algebraic  equations.  Solving  these  equations,  there  remain  ten  unknowns  (for 
(k)  (k) 

example,  ^o’  '^o  ^^i  0 »  determined  by  Eqs.  19a-b  and 

the  boundary  conditions  other  than  those  on  the  traction-free  edges. 


3.3  Complete  Solution 

Now  the  complete  solution  for  the  thermal  boundary-layer  stress  problem 
can  be  written  as 

~  1j2,4,5,6)  ,  (37a) 

u.  =  uj^^  +  u^fP^  (j  =  1,2,3)  ,  (37b) 

where  given  by  Eqs.  28a-b  and  Eqs.  31,  35, 

respectively.  The  solution  for  the  stress,  Eq.  37a,  satisfies  identically 

the  boundary  conditions,  Eqs.  12a-c,  of  the  free-edge  surface  9B  .  Along 

r 

8B  (9B  =  8B  -  3B  ) ,  residual  stresses  appear  due  to  the  particular 

S  S 


17 


solution  being  posed  in  the  previous  section.  The  residual  stress  can  be 

counterbalanced  by  the  homogeneous  solution;  thus,  on  3B  ,  one  obtains 

s 


(h) 

n 

+ 

,(h) 

n  = 

n 

,(P) 

n  , 

(38a) 

X 

X 

xy 

y 

X 

X 

xy 

y 

(h) 

n 

+ 

n  = 

n  - 

a'o) 

n  . 

(38b) 

xy 

X 

y 

y 

xy 

X 

y 

y 

(b) 

n 

+ 

,0.) 

n  = 

- 

n  - 

,(p) 

n  . 

(38c) 

xz 

X 

yz 

y 

XZ 

X 

yz 

y 

Without  orthogonality  among  eigenfunctions,  Eqs.  38a-c  may  be  satisfied 
numerically  in  a  least  square  sense  through  a  boundary  collocation  method 
by  truncating  the  eigenfunction  series.  It  is  noted  that  the  particular 
solution  is  coupled  with  the  homogeneous  one  through  the  end  conditions, 

Eqs.  19a-b.  Thus  by  matching  Eqs.  38a-c  with  the  aid  of  Eqs.  19a-b,  one  can 
determine  all  the  unknowns  explicitly.  To  this  end,  the  expressions  for  a 

z 

may  be  obtained  from  Eqs.  6d  and  37a  as 

=  -  S„,  /S„  (j  =1,2, 4, 5, 6),  (39a) 

z  3j  J  33 

=  A,  X  +  A„  y  +  A„  -  (S_.  +  a„  AT)/S__.  (39b) 

z  1  2  3  3j  3  3  33 

Complete  expressions  for  the  displacement  field  can  be  written  explicitly  in 


a  similar  manner. 
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4  .  NUMERICAL  EXAMPLES  AND  DISCUSSION 

For  illustration,  symmetric  angle-ply  composite  laminates  with  [±9]  ori- 

s 

entations  are  studied.  Thermal  loading  in  the  form  of  a  uniformly  distributed 
temperature  change  is  considered.  The  particular  laminate  configuration  and 
fiber  orientation  are  chosen  because  approximate  numerical  solutions  for  this 
problem  are  available  in  the  literature  for  comparison.  Each  individual  lamina 
is  assumed  to  be  high-modulus  graphite-epoxy  with  the  following  elastic  constants: 
Ell  =  E22  =  2.1  X  (psi).  E^2  =  20  X  10^  (psi),  G^2  =  ^13  =^23  =  0.85  x  IQ^ 

(psi),  v^2  =  ^13  =  ^93  =  0*21,  =  02  =  16  X  lo'^  (/°F),  =  0.2  x  lo“^  (/°F), 

where  the  subscripts,  1,  2  and  3,  refer  to  transverse,  thickness,  and  longitu¬ 
dinal  directions  of  the  individual  ply,  respectively.  The  and  are 
evaluated  by  using  these  material  constants,  and  the  following  relationships 
can  be  established  readily: 


=  0 

(j=l,2,3,5,6); 

=  0 

61 

(i=l,2,3,5; 

8=1,2); 

13  ij 

(l,j  i  3); 

(J=4,5,6) ; 

~(1)_~(2) 
ij  ij 

(i,j  ^4  and  i  ^  j); 

3  3 

(j=l,2); 

0  (8=1,2); 

~(1)_~(2). 
“3  “3  ’ 

(40) 

4.1  Symmetry  Cond-ittons  and  Further  Simpldoat-Lons 

The  geometric  and  lamination  symmetry  conditions  lead  to  the  following 


relationships : 

U(x,y)  =  U(x,-y),  V(x,y)  =  -V(x,-y),  W(x,y)  =  W(x,-y),  (41a-c) 

U(x,y)  =-U(-x,y),  V(x,y)  =  V(-x,y),  W(x,y)  =-W(-x,y).  (42a-c) 

Equations  41a-c  and  42a-c  may  be  written  in  equivalent  forms  as 

U  ,  (0,y)  =  V  ^(0,y)  =  W  ^  (0,y)  =  0,  (43) 

>y 

U  (x,0)  ==  V  (x,0)  =  W  (x,0)  =  0.  (44) 

^y  >y 


The  relations  provided  by  Eqs.  7  and  40  and  the  symmetry  conditions  given  in 
Eqs.  41  and  42  suggest  that 
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=  a(2>  = 
1  1 


=  0 


(i=l,2,4). 


(45) 


and  that  a  ,  a  ,  a  and  t  are  symmetric  with  respect  to  the  x  and  y  axes 
X  ^  z  xz 

and  T  and  t  are  antisymmetric  with  respect  to  the  x  and  y  axes.  Thus, 

(g) 

only  remains  to  be  determined  and  only  one-quarter  of  the  laminate  needs 

to  be  considered. 

By  using  Eqs.  40a  and  40b,  one  can  easily  show  that 


^5  -  ^5 


-  m(2) 

^3  -  m3  , 


a(^>  =  A<2) 

A3  A3 


hi.  ^  at 

:(i)  ^3  ^  :(i) 


52 


52 


(46) 

(47) 

(48) 


and  all  other  unknowns  vanish.  Thus,  the  particular  solution  takes  the 
following  forms: 


.  0, 


(49a) 


,(P)3  _ 


_51_.(1) 
;(i)  ^3 


52 


gd) 

^52 


AT, 


^(P)3  =  ^(P)3  ^  ^(P)3  ^  Q 
yz  xz  xy  ’ 


(49b) 


(49c) 


in  both  layers  (3  =  1,2).  The  displacements  may  be  shown  to  have  the  expressions, 

,(P)3  _ 


U 


AD  AD 

~(1) 

AD  53  A2 

iO) 

V 

-  12  ;(1) 

AT 

13  "  ~(1) 

^52 

;(!)  5 

®52 

X, 


(50a) 


r(p)3  _ 


AD  AD 

~(1) 

AD  53  "^22 

A^l>  + 

;(1)  _  "^22  (1) 

23  ;(i) 

A3  + 

“2  ~(1)  “5 

®52  j 

^52 

V.  / 

AT 


(50b) 
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and 

=  0.  (50c) 

Inserting  Eqs .  49a-c  into  Eq.  39b,  one  can  obtain  the  stress  component 

as 

z 


.(p)3 


1  + 


^(1) 

^32 

^53 

~(1) 

^52 

^33 

fgW)  ;(!) 

32  ^5  (1) 

~(1)  “3 

52 


AT 

,(!)• 

'33 


(51) 


Hence,  the  following  relationship  for  may  be  established: 


dx  dy  =  4 

z 


32 

^53 

~(1) 

^52 

33 

J 

b(h^  +  h^) 


+  4 


^32  “5  (1) 

~(1)  ■  “3 

^52 


b(hi  +  hp  ^ 

33 


(52) 


The  contribution  from  the  homogeneous  solution  may  be  obtained  from  Eqs.  28 
and  39a  by 


dx  dy  =  y  d  y  » 
z  ^  m  3m 

B 


m 


where 


rr 


'3m 


fo  (x,y;  6  )  dx  dy. 
3m  m 


(53) 


(54) 


Substituting  Eqs.  52  and  53  into  Eqs.  19  and  solving  the  resulting  equation 

for  A„  in  terms  of  d  give 
3  m 

=  e  +  ^  d  e  ,  (55) 

3  o  m  3m 

m 

where  ®3jjj3te  known  constants  obtained  by  integrating  associated  eigenfunctions, 
and  e^,  from  the  particular  solution. 
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The  boundary  conditions  on  (8B  -  9B  )  suggest,  by  using  Eqs .  28,  38,  43, 

r 

44,  49  and  50,  the  following  relationships  for  determining  the  constants  d  : 

m 


and 


52 


f  ^m  ^2m  ^’^’’^1’  ;(1)  S 

I  ■‘m  ^"<='••'1=  V  ■  “• 

I  (*■‘■1;  V  -  “• 

m 

J  ■‘m  "im  (''•-■'2  =  ^  ■  0. 


m 


m 


m 


I  )  -  0. 


m 


m  “2m  “2*  m 


I  ■l.  l'$P(x,-h,;S  )  -  0. 


m 


m  3  m  '  ’  2’  m- 


I  d  hf^^b,y;  6  )  =  0, 

m  im  m 

m 


m 


I '•-C  (»->'•  «»>  -o- 


m 


m 


:(i) 

’52 


AT, 


(3  =  1,2) 


(56a) 

(56b) 

(56c) 

(56d) 

(56e) 

(56f) 

(56g) 

(56h) 

(56i) 


where  h..  denote  differentiated  forms  of  g..  according  to  Eqs.  43  and  44  with 
1 J  J 

the  origin  of  '  the  Cartesian  coordinates  being  transferred  to  the  left 

free  edge  (Fig.  2).  The  boundary  conditions  shown  in  Eqs.  56a-l  are  matched 

by  a  boundary  collocation  method,  using  the  eigenfunctions  in  Eq.  28.  The 

constants  d  are  then  evaluated  by  satisfying  these  conditions  at  a  given 
m 

number  of  selected  collocation  stations.  Accuracy  and  convergence  of 
solutions  and  effects  of  collocation  points  along  the  boundaries  are 
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reported  elsewhere  [9].  For  the  numerical  results  reported  in  this  paper, 

63  terms  in  the  eigenfunction  series  and  74  collocation  stations  along  the 
boundaries  were  used  to  ensure  the  accuracy  and  convergence  of  the  solutions 
[9]. 

In  what  follows ,  laminate  thermoelasticity  solutions  determined  from 
the  current  eigenfunction  expansion  method  are  presented  first  and  compared 
with  existing  approximate  numerical  solutions  available  in  the  literature. 
Detailed  results  characterizing  the  thermal  boundary-layer  field  in  the 
symmetric  angle-ply  composites  with  various  lamination  variables  are  given 
also. 

4.2  Thermal  Boundary -Layer  Stress  Distribution 

Distributions  of  the  in-plane  and  interlaminar  thermal  stresses,  a  , 

z 

^xz’  '^y  ^yz’  the  ply  interface  of  a  [45®/-45°/-45°/45°]  graphite- 

epoxy  laminate  are  shown  in  a  semi-logarithmic  plot  in  Fig.  3.  Dotted  lines 
in  the  figure  represent  the  result  given  in  [3]  by  using  constant-strain  tri¬ 
angular  elements  (CST)  in  a  finite-element  procedure.  Solid  lines  denote  the 
present  laminate  thermoelasticity  solution  with  free-edge  stress  singularity 
being  included.  The  two  solutions  are  in  good  agreement  in  the  region  away 
from  the  laminate  boundary.  The  in-plane  stresses,  and  in  the  region 

away  from  the  edge  are  found  to  be  relatively  constant  and  to  recover  to  more 
or  less  what  classical  lamination  theory  (CLT)  predicts.  (Based  on  the 
classical  lamination  theory,  the  only  stress  induced  by  a  unit  change  of 
temperature  in  the  [±45]^  graphite-epoxy  composite  is  a  constant  in-plane 
shear  stress  =28.8  psi/°F.)  As  the  edge  is  approached,  the  difference 
becomes  gradually  appreciable  due  to  the  presence  of  the  stress  singularity, 
which  was  not  included  in  the  previous  approximate  numerical  solution.  As 
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will  be  shown  later,  within  the  boundary-layer  region  the  stress  field  is 
completely  governed  by  the  singular  terms  in  the  stress  solution,  and  the 
interlaminar  stress  increases  very  rapidly.  In  fact,  the  interlaminar  stress 
is  much  higher  than  the  in-plane  stress  component  as  the  edge  of  the 
laminate  is  Infinitesmally  approached.  Thus  thermally-induced  deformation 
and  failure  may  be  dominated  by  the  interlaminar  stresses. 


4.3  Stress  Singularity  in  Thermal  Boundary -Layer  Field 

Since  the  stress  and  strain  fields  within  the  boundary-layer  region  are 
governed  by  the  singular  terms  in  the  present  laminate  thermoelasticity  solu¬ 
tion,  theoretically,  the  thermally  induced  stresses  determined  from  the  solution 
are  unbounded  at  the  intersection  of  the  ply  interface  and  the  laminate  edge. 
Thus  the  near-field  stress  may  be  expressed  in  a  general  form  as 


a. 

1 


3 

I 


k=l 


+  D 


i(k+3)^k 


+  0 (higher-order,  non-singular  terms) 


(i=l,2,3 . 6)  ,  (57) 

where  Z  and  Z  have  their  origin  at  the  intersection  of  the  ply  interface  and 
the  edge  of  the  laminate  (Fig.  2);  6^^  is  the  order  of  the  thermal  boundary- 
layer  stress  singularity  which  is  the  smallest  eigenvalue  satisfying 


0  >  Re[6^]  >  -1  (58) 

among  all  the  5^  determined  from  the  characteristic  equation  [9]  in  solving 
for  the  homogeneous  solution  of  the  governing  partial  differential  equations. 
The  order  of  the  boundary-layer  stress  singularity  is  noted  [9]  to  depend 
only  upon  lamina  constitutive  properties  and  fiber  orientations  of  the  adjacent 
plies.  Numerical  results  of  the  first  four  nonzero  eigenvalues  6^  for  sym¬ 
metric  angle-ply  [0/-0/-0/0]  graphite-epoxy  laminates  are  given  in  Table  1. 
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Note  that  6^  corresponds  to  the  order  of  boundary-layer  stress  singularity 
and  that  zero  and  intergers,  n,  are  always  eigenvalues  for  the  problem.  In 
this  composite  system,  6^  has  a  value  which  is  rather  weak  as  compared  with 
other  typical  singular  stress  problems  such  as  an  elastic  crack  problem. 

Among  various  6  studied,  it  appears  that  the  [±51“]^  graphite-epoxy  laminate 
has  the  strongest  stress  singularity  [9]. 

4.4  Thermal  BoundaTy -Layer  Stress  Intensity  Factors 

For  a  composite  laminate  with  given  fiber  orientations,  the  coefficients 
of  the  singular  terms  in  Eq.  57  characterize  amplitudes  of  the  thermal  stress 
and  strain  in  the  boundary-layer  region.  Because  the  boundary-layer  stresses 
are  most  crucial  along  the  ply  interface,  i.e.,  the  x-axis,  and  become  singular 
at  the  Interface/edge  Intersection,  it  is  possible  to  define  the  amplitudes 
of  the  singular  thermal  boundary-layer  stresses  by 

—6 1 

K.  =  lim  X  (i=l,2,3 . 6).  (59) 

^  x->0  ^ 

The  are  dependent  upon  geometric  variables  of  the  composite  (e.g.,  ply 
thickness,  number  of  layers),  lamination  parameters  (e.g.,  fiber  orientation, 
stacking  sequence),  and  mechanical  and  thermal  loading  conditions.  The 
fundamental  structure  of  the  thermal  boundary-layer  stress  solution  shown 
in  Eqs.  57  and  59  resembles  that  of  an  elastic  crack  problem  (except  that  6^ 
has  a  value  of  -0.5  in  the  crack-tip  stress  field).  Also,  the  nature  of 
is  similar  to  the  so-called  crack-tip  stress  intensity  factors  in  linear 
elastic  fracture  mechanics.  Thus,  in  this  context,  it  may  be  appropriate  to 
denote  as  "thermal  boundary-layer  stress  intensity  factors"  or  "thermal 
free-edge  stress  intensity  factors"  for  the  composite  laminate.  Values  of 
for  the  [0/-0/-e/0]  graphite-epoxy  composite  with  all  laminae  being  of 
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equal  thicknesses  under  unit  thermal  loading  are  determined  in  Table  2,  in 

-<Si 

which  carry  the  unit  of  [psi-in.] 

The  associated  with  the  interlaminar  stresses  are  found,  in  general, 
significantly  larger  than  those  associated  with  the  in-plane  stress  components. 


The  dominance  of  the  interlaminar  thermal  stresses  a  and  x  in  the  boundary- 

y  yz 

layer  region  shown  in  Fig.  3  is  clearly  illustrated  by  the  high  values  of  K2 
and  K^,  which  are,  in  fact,  one  or  two  orders  of  magnitude  higher  than  the 
associated  with  a^,  0^  and  x^^  for  all  6  studied.  The  high  negative  value  of 


K2  indicates  that  a  large  compressive  interlaminar  normal  stress  is 
developed  near  the  edge.  Note  that  is  found  to  vanish  for  all  angle-ply 
[±9]^  composite  laminates,  due  to  the  symmetry  of  ply  orientations  and  traction- 


free  edge  conditions. 


4.5  Through-Thickness  Distribution  of  Thermal  Boundary-Layer  Stresses 

The  unique  features  of  the  thermal  boundary-layer  effect  are  further 

illustrated  by  through-the-thickness  distributions  of  in-plane  and  interlaminar 

thermal  stresses  near  the  laminate  boundary.  The  in-plane  thermal  stress 

in  the  thickness  direction  at  different  distances  away  from  the  edge  are  shown 

in  Fig.  4.  The  current  solution  is  in  agreement  with  previous  results  [3] 

that  0^  is  compressive  in  the  most  part  of  the  section  near  the  laminate 

boundary  except  for  the  region  closest  to  the  ply  interface  where  boundary- 

layer  stress  singularity  dominates.  Through-the-thickness  distributions  of 

the  most  dominant  thermal  interlaminar  stress  x  are  given  in  Fig.  5.  The 

yz 

gradient  of  x^^  in  the  y-directlon  increases  rapidly  as  the  laminate  edge  is 
approached.  Again  the  laminate  elasticity  solution  in  the  boundary-layer 
region  differs  from  the  approximate  solution  [3]  near  the  interface  but  in 
in  good  agreement  in  the  far  field.  The  next  two  figures  (Figs.  6  and  7) 
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provide  information  on  through- the- thickness  distributions  of  thermal  inter¬ 
laminar  shear  and  normal  stresses,  and  in  the  [45°/-45°/-45°/45°] 
graphite-epoxy  laminate.  The  is  found  to  change  significantly  with  the 

thickness  coordinate,  alters  its  sign  across  the  interface,  and  reaches  a 
maximum  value  at  a  small  distance  y  away  from  the  interface.  Within  the 
boundary-layer  region,  becomes  very  small  along  the  interface  and  vanishes 
identically  at  the  laminate  boundary  as  envidenced  by  =  0.  The  distribution 
of  0  at  several  distances  near  the  boundary  of  the  laminate  is  shown  in 

y 

Fig.  7.  At  any  given  x/b,  is  generally  very  small  at  a  distance  y  away 
from  the  interface,  and  reaches  a  higher  level  as  the  interface  y  =  h  is 
approached.  As  one  moves  towards  the  edge,  the  interlaminar  normal  stress 
becomes  very  significant  in  compression.  As  x  0  and  y  h  the  stress 
solution  becomes  unbounded  due  to  the  stress  singularity  at  the  intersection 
point. 

4.6  Thermal  Boundary -Layer  Width 

The  rapid  increase  of  thermal  stress  has  been  observed  to  be  restricted 
to  within  a  very  localized  region  near  the  edge  of  the  laminate — the  so-called 
"thermal  boundary-layer  width."  The  thermal  stresses  developed  in  the 
boundary- layer  region  are  inherently  three-dimensional  and  cannot  be  determined 
by  classical  lamination  theory.  The  singular  nature  and  the  extent  of  pertur¬ 
bation  of  the  thermal  boundary  stress  are  considered  to  be  of  vital  importance 
in  controlling  initiation  of  interlaminar  fracture  (or  delamination)  and 
strength  degradation.  The  extent  of  perturbation  of  the  laminate  thermal 
stress  field  can  be  characterized  by  the  boundary-layer  width  (or  thickness),  B. 
Pipes,  et  al.  [10]  defined  the  boundary-layer  thickness  as  the  distance 
from  the  edge,  at  which  the  interlaminar  stress  is  about  3  per  cent  of 
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the  value  calculated  at  the  intersection  of  the  ply  interface  and  the  edge  of 
the  laminate.  The  validity  of  this  definition  is  somewhat  questionable 
because  the  interlaminar  stresses  are  singular  at  this  point.  In  this 
study  an  alteimative  definition  of  the  boundary-layer  width  is  proposed  on 
the  basis  of  the  strain  energy  density  distribution  in  the  composite  laminate. 

The  strain  energy  density  distributions  E(x,y)  along  the  interface  of 
the  [0/-0/-9/6]  graphite-epoxy  laminates  are  shown  in  Fig.  8.  The  strain 
energy  density  remains  relatively  constant  in  the  far  field  where  classical 
lamination  theory  holds,  and  increases  drastically  by  an  order  of  magnitude 
as  the  edge  is  approached.  In  this  paper,  the  boundary-layer  width  B  in  a 
composite  laminate  is  defined  as  the  distance  away  from  the  edge  where  the 
strain  energy  density  along  the  interface  is  three  per  cent  higher  than  the 
nominal  value  obtained  in  the  far  field.  In  general,  E(B,h'*’)  differs 
slightly  from  E(B,h  )  due  to  the  discontinuous  in-plane  stress  components 
at  y  =  h"*"  and  h  ;  thus,  an  average  value  of  B  is  designated  as  the  width  of 
the  boundary-layer  region.  Based  on  this  definition,  values  of  B/W  for 
[0/-0/-0/0]  graphite-epoxy  laminates  are  evaluated  and  shown  in  Fig.  9 
It  is  obvious  from  the  figure  that  the  [45'’/-45°/-45‘’/45'’]  graphite-epoxy 
laminate  has  a  higher  value  of  B/W  than  those  composites  with  other  fiber 
orientations.  As  0  changes  towards  either  direction,  B/W  decreases  rapidly. 
When  0  has  an  angle  of  [±0°]  or  [±90° ],  B/W  vanishes  identically  indicating 
that  there  is  no  boundary- layer  effect  in  these  cases  since  the  two  adjacent 
plies  are  identical. 

4.7  Effects  of  Ftbev  Ortentation 

Effects  of  fiber  orientation  on  the  thermal  boundary-layer  response  in 
[0/-0/-0/0]  graphite-epoxy  laminates  are  best  illustrated  by  boundary-layer 
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thermal  stresses  along  0/-0  ply  interface.  The  thermal  stresses  are  observed 
to  be  significantly  affected  by  the  alternation  of  fiber  orientation.  Figures 
10  and  11  depict  the  distributions  of  thermal  interlaminar  shear  stresses,  t 

xy 

and  y  =  h.  The  Is  relatively  small  for  all  fiber  orientations 

studied,  and  the  change  of  fiber  orientation  only  alters  the  amplitude  of  the 

Interlaminar  shear  stress  slightly.  The  x  reaches  its  maximum  before  it 

xy 

vanishes  at  the  laminate  boundary,  where  traction-free  boundary  conditions 
are  satisfied  exactly.  The  other  thermal  interlaminar  shear  component  x^^  is 
more  significant  in  the  boundary-layer  region  than  The  x^^  has  a  higher 

value  along  the  interface  of  the  [45°/-45°/-45°/45°]  graphite-epoxy  laminate 
than  those  in  other  ply  ocnf iguratlons  due  to  the  higher  value  of  6^.  As  is 
expected,  the  interlaminar  shear  stress  becomes  unbounded  as  the  laminate 
boundary  is  approached  due  to  the  stress  singularity.  As  0  moves  towards 
either  side,  the  amplitude  of  the  interlaminar  shear  stress  decreases.  The 
distribution  of  interlaminar  normal  stress  along  the  interface  of  angle-ply 
laminates  with  different  0  is  shown  in  Fig.  12.  The  a  is  found  to  be 

y 

vanishingly  small  in  the  far  field.  It  is  small  in  tension  first,  then 
changes  its  sign,  and  becomes  compressive  as  the  free  edge  is  approached. 

Like  the  interlaminar  shear  stress,  becomes  unbounded  at  the  free  edge  (i*e., 
at  X  =  b  and  y  =  h)  and  is  significantly  only  within  the  boundary-layer  region. 
In  Fig.  13,  distributions  of  In-plane  thermal  stress  in  [0/-0/-0/0] 
graphite-epoxy  are  given.  The  laminate  thermoelasticity  solution  reveals 
that  has  a  very  small  value  in  the  far  field,  and  remains  relatively 
constant  before  rapid  increase  in  its  magnitude,  as  the  laminate  boundary 
is  approached.  This  is  apparently  different  from  the  prediction  of  classical 
lamination  theory,  which  suggests  that,  in  the  symmetric  angle-ply 
[0/-0/-0/0]  composite  laminate,  the  in-plane  thermal  stress  vanishes 
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throughout  the  composite.  At  the  Intersection  of  the  free  edge  and  ply 
interface,  again  has  an  unbounded  value  due  to  the  stress  singularity 
at  that  point.  In  the  boundary-layer  region,  along  the  interface  in  the 
[±45]^  laminate  is  found  to  have  a  higher  value  than  in  the  composites  with 
other  ply  configurations.  As  0  changes  from  45°,  the  magnitude  of  de¬ 
creases  appreciably.  In  fact,  the  results  shown  in  Fig.  13  indicate  that  0^ 
becomes  negligibly  small  for  the  cases  with  6  >  75°. 

4.8  Effects  of  Relative  Ely  Thickness 

Another  important  laminate  variable  investigated  in  this  study  is  the 
effect  of  ply  thickness  or  volume  of  the  layer  that  is  stressed  inter- 
lamlnarly,  since  transverse  deformation  and  failure  in  composites  are  affected 
significantly  by  lateral  constraints  in  the  laminates.  It  has  been  found 
that  deformation  and  fracture  in  certain  composite  systems  change  significant¬ 
ly  with  the  ply  thickness  while  stacking  sequence  remains  the  same.  In 
this  paper,'  influences  of  ply  thickness  on  the  thermal  boundary-layer  stress 
in  [45°/-45°/-45°/45°]  graphite-epoxy  laminate  with  various  h^^/W  are  examined 
(W  being  kept  constant) .  Following  the  aforementioned  analytical  procedure 
and  solution  scheme,  a  parametric  study  on  thermal  boundary-layer  response  in 
the  composite  laminates  with  various  ply  thicknesses  h^/W  has  been  conducted. 
Numerical  results  showing  the  effect  of  h^/W  on  thermal  boundary-layer 
stress  intensity  factors  are  given  in  Table  3.  The  results  indicate  that, 
for  a  given  laminate  configuration  and  fiber  orientation,  the  change  of  h^/W 
does  not  alter  the  thermal  stress  singularity  but  does  affect  thermal  boundary- 
layer  stress  intensity  factors  appreciably.  The  composite  laminate  with 
h^/W  -  0.4  ^  0.5,  seems  to  have  smaller  values  of  than  other  cases  studied. 
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As  h^/W  approaches  0  or  1,  the  boundary-layer  stress  intensity  factors  reach 
their  maxima.  The  higher  values  of  developed  in  the  composite  by  changing 
ply  thickness  h^/W  provide  an  important  basis  for  evaluation  of  initiation  of 
transverse  cracking  and  interply  delamination  in  composite  laminates. 
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5.  SUMMARY  AND  CONCLUSIONS 

A  study  of  thermal  boundary-layer  stresses  in  composite  laminates  has 
been  presented.  Formulation  of  the  problem  is  based  on  the  theory  of 
anisotropic  laminate  thermoelasticity.  With  the  aid  of  Lekhnitskii ' s 
comp lex- variable  stress  functions,  an  eigenfunction  expansion  method  is 
used  to  establish  a  system  of  coupled,  governing  partial  differential 
equations  for  the  problem.  Numerical  results  for  S3nnmetric  angle-ply 
[0/-6/-6/0]  graphite-epoxy  laminates  are  obtained.  Effects  of  lamination 
and  geometric  variables  on  the  thermal  boundary-layer  stress  singularity 
and  distributions  are  studied.  Based  on  the  information  discussed  in 
the  previous  sections,  the  following  conclusions  may  be  reached: 

1.  Thermal  stresses  in  the  boundary-layer  region  of  a  composite 
laminate  are  inherently  three  dimensional  in  nature.  They 
cannot  be  calculated  by  the  classical  lamination  theory,  but 
can  be  determined  explicitly  by  the  current  approach. 

2.  The  thermal  stress  field  in  the  boundary-layer  region  is 
singular  in  general.  By  using  an  eigenfunction  expansion 
method,  one  can  determine  the  order  of  the  boundary-layer 
stress  singularity  by  solving  the  characteristic  transcendental 
equation  obtained  from  the  homogeneous  solution.  The  order 

of  the  boundary-layer  stress  singularity  depends  on  anisotropic 
thermoelastic  properties  of  adjacent  plies  in  the  composite. 

3.  The  boundary-layer  thermal  stress  field  may  be  characterized  by 
"thermal  boundary-layer  stress  intensity  factors"  or  "thermal 
free-edge  stress  Intensity  factors."  The  are  functions  of 
anisotropic  thermoelastic  constants  of  laminae,  ply  orientation 
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and  laminate  geometry,  and  may  be  used  to  evaluate  strength 
degradation  and  initiation  of  interlaminar  fracture 
(delamination)  and  transverse  cracking.  Their  values  can 
be  determined  by  various  methods  such  as  the  boundary 
collocation  method  and  the  finite  element  method. 

4.  Thermal  boundary-layer  thickness  which  characterizes  the  domain 
where  classical  lamination  theory  does  not  hold  can  be 
determined  explicitly  by  considering  the  change  of  strain 
energy  density  along  the  ply  interface.  The  thermal  boundary- 
layer  thickness  depends  on  lamination  variables ,  geometric 
parameters,  thermal  loading  conditions  and  thermoelastic  ply 
properties.  In  l-0]g  graphite-epoxy  composites,  the  case  of 

9  =  45°  possesses  a  higher  thermal  boundary-layer  thickness 
than  those  of  other  fiber  orientations  studied. 

5.  In  comparison  with  previous  approximate  solutions,  good 
agreement  in  the  far  field  is  observed,  but  appreciable 
discrepancy  near  the  laminate  boundary  occurs.  iThe  difference 
is  attributed  to  the  fact  that,  in  the  boundary-layer  region, 
the  thermoelastic  solution  is  completely  governed  by  the 
singular  terms  which  previous  approximate  solutions  failed 

to  include. 

6.  The  current  method  of  approach  is  also  valid  for  asymmetric 
composite  laminates,  since  the  bending,  twisting  and  rotational 
components  of  deformation  are  included  in  the  formulation. 
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Table  1 

First  four  nonzero  eigenvalues  for  thermal  boundary-layer  stresses 
in  sjnnmetric  angle-ply  [6/-0/-6/6]  graphite-epoxy  laminates* 


0“ 

^1 

6 

2,3 

^4 

15° 

-0.64322 

E-3 

0.99670  ± 

0.04191  i 

1. 

o 

o 

-0.11658 

E-1 

0.95521  ± 

0.15271  i 

1. 

45° 

-0.25575 

E-1 

0.88147  ± 

0.23401  i 

1. 

o^ 

o 

o 

-0.23346 

E-1 

0.83074  ± 

0.27138  i 

1. 

75° 

-0.89444 

E-2 

0.86469  ± 

0.25007  i 

1. 

*6-  corresponding  to  the 

strength 

of  thermal 

boundary-layer 

stress 

singularity. 
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Table  2 

Thermal  boundary-layer  stress  intensity  factors,  K^,  along  the 
interface  of  [0/-0/-6/0]  graphite-epoxy  composite  laminates*”^ 

e  K3  K3 

15°  4.5996  E-1  -3.8534  E  2  -3.7120  E  1  -5.2624  E  3  8.0636  E  0  0 

30°  4.7031  E  0  -1.6466  E  2  8.6078  E  0  -5.7404  E  2  2.2357  El  0 

45°  1.4546  E  1  -1.9130  E  2  -7.4958  E-1  -4.1747  E  2  3.6662  El  0 

60°  1.3746  E  1  -1.7881  E  2  -2.2426  E  1  -3.7858  E  2  2.4394  El  0 

75°  4.5752  E  0  -1.5101  E  2  -2.9205  E  1  -4.9611  E  2  8.5717  E  0  0 

*h3  =  h^  =  h,  b  =  8h 

"^Values  of  are  per  °F  change 
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Table  3 

Lamina  thickness/ volume  effects  on  thermal  boundary-layer  stress 
intensity  factors  for  [45“/-45°/-45°/45°]  Graphite -Epoxy  Laminates*^ 


h^/W 


0.2  1.4870  E  1  -1.9554  E  2  -7 
0.4  1.4509  E  1  -1.9079  E  2  -7 
0.5  1.4548  E  1  -1.9130  E  2  -7 
0.6  1.4707  E  1  -1.9339  E  2  -7 
0.8  1.5580  E  1  -2.0487  E  2  -8 


6620 

E-1 

-4.2672 

E 

2 

3.7475 

E 

1 

0 

4759 

E-1 

-4.1636 

E 

2 

3.6565 

E 

1 

0 

4958 

E-1 

-4.1747 

E 

2 

3.6662 

E 

1 

0 

5780 

E-1 

-4.2204 

E 

2 

3.7064 

E 

1 

0 

0275 

E-1 

-4.4708 

E 

2 

3.9263 

E 

1 

0 

*2b  =  8W  where  W  =  constant  =  half  laminate  thickness 
"^Values  of  are  per  “F  change 
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Fig.  3  In-plane  and  Interlaminar  Thermal  Stresses  along  Interface  in 

[45*/-45V-45®/45«]  GrAPHITE-EpOXY  LaMINATE^  (h,  =h2=h,  b  =  8h) 
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ct^/AT  (psi/^F) 


Fig.  ^  Through-Thickness  Distribution  of  Boundary-Layer  Thermal 
Stress  in  [45V-45V-45V45®]  Graphite-Epoxy  Composite 

(h^  =h2  =h,  b=8h) 
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Fig.  5  Through-Thickness  Distribution  of  Boundary-Layer  Thermal 
Stress  Ty^  in  [45‘’/-45V-45“/45°]  Graphite-Epoxy  Composite 
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Fig.  6  Through-Thickness  Distribution  of  Boundary-Layer 
Thermal  Stress  x^y  in  [45V-45V-45V45‘’]  Graphite- 
Epoxy  Composite 
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Fig.  7  Through-Thickness  Distribution  of  Boundary-Layer 
Thermal  Stress  in  [45V-45V-45V45°]  Graphite- 
Epoxy  Composite 
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/AT  (psi/^F) 


Fig.  10  Therm. 
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Fig.  li  Thermal  Interlaminar  Shear  Stress  Ty^  along  Inter¬ 
face  IN  [e/-0/-e/e]  Graphite-Epoxy  Composites 
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Fig.  12  Thermal  Interlaminar  Normal  Stress  along  Inter¬ 
face  IN  [e/-e/-0/e]  Graphite-Epoxy  Composites 
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